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ABSTRACT 

We study the temporal and coarse spectral properties of 268 bursts from SGR 
1806-20 and 679 bursts from SGR 1900+14, all observed with the Rossi X-Ray 
Timing Explorer/Proportional Counter Array. Hardness ratios and temporal 
parameters, such as Tgo durations and rgo emission times are determined for 
these bursts. We find a lognormal distribution of burst durations, ranging over 
more than two orders of magnitude: Tgo ~ 10^^ to > 1 s, with a peak at ~ 0.1 
s. The burst light curves tend to be asymmetrical, with more than half of all 
events showing rise times tr < 0.3 Tgo. We find that there exists a correlation 
between the duration and fiuence of bursts from both sources. We also find a 
significant anti-correlation between hardness ratio and fiuence for SGR 1806—20 
bursts and a marginal anti-correlation for SGR 1900+14 events. Finally, we 
discuss possible physical implications of these results within the framework of 
the magnetar model. 

Subject headings: X-rays: bursts - gamma rays: bursts - stars: individual (SGR 
1806-20) - stars: individual (SGR 1900+14) 



1. Introduction 

Soft gamma repeaters (SGRs) are a small class of objects that are characterized by 
brief and very intense bursts of soft gamma-rays and hard X-rays. They are distinguished 
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from classical gamma-ray bursts (GRBs) by their repeated periods intense activity, during 
which dozens of bursts with energies approaching 10^^ ergs are recorded. SGR bursts have 
significantly softer spectra than classical GRBs; the former are being well fit by an optically- 
thin thermal bremsstrahlung model with temperatures kT = 20 — 40 keV. Two SGRs (0526- 
66 and 1900+14) have emitted one giant flare each: events that are much more energetic 
[E ~ ]^q44-45 gj,g^ contain a very hard spectral component within the first ~ 1 s (Mazets 
et al. 1979; Cline et al. 1981; Hurley et al. 1999; Mazets et al. 1999; Feroci et al. 1999). 
For a review of the burst and persistent emission properties of SGRs, see Hurley (2000). 

In 1992, it was suggested that SGRs arc strongly magnetized {B > 10^^ G) neutron stars, 
or magnctars (Duncan & Thompson 1992, see also Paczynski 1992). This model suggests 
crustquakes as a plausible trigger for the short SGR bursts (as well as the giant flares): a 
sudden fracture of the rigid neutron star crust, driven by the build-up of crustal stresses as 
the strong magnetic fleld gradually diffuses through the dense stellar matter (Thompson & 
Duncan 1995 [hereafter TD95]; 1996). The motion of the crust shears and twists the external 
magnetic fleld, and in the process releases both elastic and magnetic energy. 

Cheng et al. (1996) studied a set of 111 SGR 1806-20 bursts (detected with the In- 
ternational Gometary Explorer, IGE) and determined that some of their properties, such as 
size and cumulative waiting-time distributions are similar to those of earthquakes. Recently, 
Gogiis et al. (2000) confirmed these similarities for SGR 1806—20 using a larger sample of 
401 bursts (111 detected with the Burst and Transient Source Experiment, BATSE, aboard 
the Compton Gamma Ray Observatory, CGRO; 290 with the the Proportional Counter Ar- 
ray, PGA, on the Rossi X-ray Timing Explorer, RXTE). The same similarities were found for 
SGR 1900+14 when we analyzed a sample of 1024 bursts (187 observed with CGRO/BATSE; 
837 with the RXTE/PCA) (Gogii§ et al. 1999). Furthermore, we reported evidence for a cor- 
relation between duration and fluence for the SGR 1900+14 bursts, similar to the correlation 
seen between the duration of strong ground motion at short distances from an earthquake 
region and the seismic moment (oc energy) of an earthquake (Lay & Wallace 1995). This 
similarity strongly suggests that SGR bursts, like earthquakes, may be manifestations of self- 
organized critical systems (Bak, Tang & Wiesenfeld 1988), lending support to the hypothesis 
that SGR bursts involve the release of some form of stored potential energy. 

Statistical studies such as the above provide important clues not only about the mecha- 
nism by which energy is injected during an SGR burst, but also the mechanism by which it is 
radiated. For example, in the magnetar model the stored potential energy may be predom- 
inantly magnetic, and the electrical currents induced by crustal fractures may be strongly 
localized over small patches of the star's surface, somewhat as they are in solar flares. In 
such a situation, the energy density in the non-potential magnetic fleld is high enough that 
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excitation of high frequency Alfven motions leads to rapid damping and the creation of a 
trapped fireball (through a cascade to high wavenumber: Thompson & Blaes 1998). By 
contrast, if energy is injected more globally, then an approximate balance between injection 
and radiation will occur in the bursts of luminosity less than ~ 10^^ erg s~^ (Thompson et 
al. 2000). 

Here, wc investigate detailed temporal characteristics of a large subset of SGR 1806—20 
and SGR 1900+14 bursts observed with the RXTE/PGA during their burst active periods 
in 1996 and 1998, respectively. We apply to the SGR bursts some temporal analysis methods 
that were originally developed for the study of GRBs. We also study the spectral variations 
of the bursts as a function of the burst fluence and duration. In § 2, we briefly review the 
RXTE/PCA observations. Section 3 describes the data analysis techniques and our results 
are discussed in 5 4. 



2. Observations 

The PGA instrument (Jahoda ct al. 1996) consists of five Xe proportional counter units 
(PGUs) sensitive to energies between 2—60 keV with 18% energy resolution at 6 keV. The 
PGA has a total effective area of ~ 6700 cm^. 

SGR 1806-20 : The RXTE/PGA observations of SGR 1806-20 were performed 
during a burst active period of the source in 1996 (between November 5 and 18) for a total 
effective exposure time of 136.8 ks. Using the burst search procedure described in Gogii§ 
et al. (2000), we have identified 290 bursts from the source. The number of integrated 
counts (2—60 keV) for these bursts range between 22 and 34550. Using the count-to-energy 
conversion factor given in Gogii§ et al. (2000) this corresponds to burst fiuences between 
1.2 X 10~^'^ and 1.9 x 10~^ ergs cm~^ (E > 25 keV). Assuming isotropic burst emission, the 
corresponding energy range is 3.0 x 10^^ - 4.9 x 10^^ ergs, for a distance to SGR 1806—20 of 
14.5 kpc (Gorbel et al. 1997). To facilitate our analysis, we chose events that were clustered 
together during two very active epochs of the source, resulting in 268 events recorded on 5 
(MJD 50392) and 18 (MJD 50405) November 1996. 

SGR 1900+14 : RXTE observations of SGR 1900+14 took place between 1998 June 
2 and December 21, for a total effective exposure time of 224.1 ks. Using the same burst 
search algorithm as before, we identified a total of 837 bursts from the source with integrated 
counts ranging between 22 and 60550. The fiuences of these bursts range from 1.2 x 10"^° 
to 3.3 X 10"''' ergs cm"^ (E > 25 keV) corresponding to an energy range of 7 x 10^^ - 2 x 10^^ 
ergs (assuming isotropic emission at 7 kpc [Vasisht et al. 1994]). Similarly, we selected 679 
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events which occurred during a very active period of the source between 1998 August 29 
(MJD 51054) and September 2 (MJD 51058). 

3. Data Analysis and Results 

3.1. Duration (Tgo) estimates 

Originally defined for cosmic GRBs, the Tgo duration of a burst is the time during 
which 90% of the total (background-subtracted) burst counts have been accumulated since 
the burst trigger (Kouveliotou et al. 1993). We calculated the Tgo duration of SGR bursts 
using event-mode PGA data (2-60 keV) with 1/1024 s time resolution. For each burst we 
collected the cumulative counts for an 8 s continuous stretch of data starting 4 s before its 
peak, tp. We then fit the cumulative count distribution between two user-selected background 
intervals ^ with a first order polynomial plus a step function (for the burst) and subtracted 
the background counts. By using a first order polynomial to fit the cumulative counts, we 
assume the PGA background remains flat over each 8 s segment. The resulting height of 
the step function gives the total burst counts. Figure 1 depicts the steps of the Tgo estimate 
procedure for one of the SGR 1900-1-14 bursts. 

As the selected bursts occurred during extremely active periods (for both sources), 
there were quite a number of cases where bursts were clustered very close together. During 
these active episodes, many crustal sites on the neutron star may be active, releasing stored 
potential energy in the form of bursts with a large variety of time profiles. Hence, it is 
important to distinguish single pulse events from events with multiple peaks (which may 
involve multiple fracture sites). To do this, we applied arbritrary but consistent criteria to 
our data. We classified an event as multi-peaked if the count rate at any local minimum 
(in 7 ms time bins) is less than half the maximum value attained subsequently in the burst 
(see Figure 2, middle plots). Otherwise the event was classified as single-peaked (Figure 2, 
top plots). When the count rate dropped to the noise level between peaks, the event was 
classified as a single, multi-peaked burst if and only if the time between peaks was less than 
a quarter of the neutron star rotation period ( 1.3 s for SGR 1900+14, 1.9 s for SGR 1806-20; 
see Figure 2, bottom plots). We found that 262 of the 679 bursts from SGR 1900+14, and 
113 of the 268 SGR 1806-20 bursts were multi-peaked by these criteria. Note that the total 
counts per burst (or " count fluences" ) were found to span nearly the same range in the two 



^Approximately equal background intervals were selected before and after the burst, except in cases where 
separation between bursts was too small; there we chose smaller, but significant, post/pre burst intervals. 



- 5 - 



classes of bursts, in both sources. This suggests that the hght curve peak structure does not 
correlate strongly with total energy. 

In the context of cosmic GRBs, the value of T90 can be systematically underestimated in 
faint events due to low signal-to- noise (Koshut et al. 1996). We investigated this effect in our 
SGR analysis using extensive numerical simulations. We created three time profiles based 
upon the observed time profiles; a single two-sided Gaussian (with right-width twice the left 
width), two two-sided Gaussians whose peaks are separated by 0.5 s and the peak rate of 
the second pulse is 0.6 of the first one, and two two-sided Gaussians with peak separation 
of 1 s. For each profile, we varied the peak rates (eight vahies between 2600 counts s~^ 
and 55000 counts s"^) and the width of the Gaussians (eight values from 12 ms to 120 
ms). We determined the T90 duration for each combination before adding noise. Then for 
each combination, we generated 800 realizations including Poisson noise and determined the 
respective T90 values for these simulations. We found that the fractional difference, FD = 1 
— (T9o,s / Tgo^a) bctween the actual value, Tgo^a, (in the absence of noise) and the simulated 
one, T9o,s, has a strong dependence on the peak rate (i.e. signal-to-noise ratio), such that 
FD can be as high as 0.22 at count rates of 2600 counts s^^. It is, however, less than 0.05 
for count rates greater than 4800 counts s~^. The variation of FD with respect to the pulse 
profile is insignificant. Using the results of our simulations we obtained a T90 correction as 
a function of peak rate. We estimated the peak rate of each burst using a box-car averaging 
technique with a box width of 1/512 s. We then corrected T90 values for the S/N effect using 
our correction function and these peak rates. 

Our final data set thus comprises 455 SGR 1900+14 bursts with corrected T90 dura- 
tions between 9 ms and 2.36 ms (Figure 3, solid histogram). We were unable to determine 
statistically significant Tgo durations for 187 SGR 1900-1-14 bursts because there were less 
than ~ 40 counts per event. The sohd curve in Figure 3 is a best fit log-Gaussian function 
to all T90 values which peaks at 93.9 ± 0.2 ms {a — 0.35 ± 0.01; where a is the width of the 
distribution in decades). The dashed histogram in Figure 3 is the distribution of T90 values 
of single pulse bursts whose log-Gaussian mean is 46.7 ± 0.1 ms {a = 0.21 ± 0.01). Also in 
Figure 3, the dash-dot histogram displays the T90 distribution of multi-peaked bursts and 
it peaks at 148.9 ± 0.2 ms {a = 0.26 ± 0.02). For SGR 1806-20 bursts, we determined 
corrected T90 values of 190 bursts which range between 16 ms and 1.82 s (49 bursts were too 
weak to determine their T90 values). The distribution of all T90 values is shown in Figure 
4 (sohd histogram). A log-Gaussian fit to this distribution yields a peak at 161.8 ± 0.2 ms 
{a — 0.34 ± 0.02). Similarly, a log-Gaussian fit to the T90 distribution of the single pulse 
bursts (dashed lines in Figure 4) peaks at 88.1 ± 0.1 ms {a = 0.19 ± 0.03) and a fit to the 
multi-peaked bursts (dash-dot lines in Figure 4) yields a peak at 229.9 ± 0.3 ms {a — 0.32 
± 0.03). 
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In order to quantify the time profile symmetry of SGR bursts we determined the ratio 
of their rise times (t^, i.e. the interval between T90 start time and peak time tp) to T90 
durations of the single pulse events. In Figure 5, we plot the distributions of these ratios 
for bursts with single pulse (solid histogram) and milti-peaks (dashed histogram) from SGR 
1900+14 (left) and those from SGR 1806—20 (right). The majority of single pulse events 
from both sources have t^ / T90 values less than 0.3, showing that SGR bursts decay slower 
than they rise, or in other words that SGR pulse profiles are asymmetric. We find that the 
average values os the t^ / T90 ratios are 0.29 and 0.27 for SGR 1900+14 and SGR 1806-20, 
respectively. Remarkably coincident average values, along with very similar distributions of 
tr I T90 ratios of both sources suggest a similarity of the asymmetry in the temporal profiles 
of SGR bursts. 



3.2. Emission time (r9o) estimates 

Emission time, tjv, was introduced for cosmic GRBs as a complimentary temporal 
parameter to T90 (Mitrofanov et al. 1999). Emission time is the time over which a fixed 
percentage, N%, of the total burst emission was recorded, starting from the peak of the 
event and moving downward in flux. For each burst, we determined the average background 
level from the background intervals used in the T90 estimating procedure. The background- 
subtracted count bins within the burst interval (i.e. from the end of prc-burst background 
range to the beginning of post-burst background range) were then ordered by decreasing 
count rate. Starting with the highest count rate bin (i.e. the peak), we added the counts of 
each successively weaker bin until 90% of the total burst counts were accumulated. The rgo 
emission time of the burst is then the total time spanned by the accumulated bins. 

As with the T90 parameter, when the S/N ratio of the burst is low, the measured emission 
time can be systematically smaller than the actual value (Mitrofanov et al. 1999). In order 
to correct for this systematic error, we constructed numerous simulated profiles identical to 
those described in the previous section. We found that, similar to T90 estimates, the FD 
between the actual and simulated rgo emission times is strongly dependent on the peak rate. 
FD is ~ 0.26 at count rates of 2600 counts s"-*^ and less than 0.06 for count rates greater 
than 4800 counts s~^. Similar to the T90 measurement, this effect is only weakly dependent 
on the temporal profile of the burst. Using our rgo correction function and the peak rates, 
wc obtained the corrected rgo emission times. The corrected values range between 4.6 and 
412.8 ms for SGR 1900+14 bursts, and between 4.8 and 559.2 ms for SGR 1806-20 bursts. 
In Figure 6, we show the distributions of rgo emission times for SGR 1900+14 (/e/t) and 
SGR 1806—20 (right). The dashed lines in both plots are the best fit log-Gaussian curves 
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which peak at 49.6 ± 0.1 ms {a = 0.28 ± 0.01) for SGR 1900+14 events, and at 82.3 ± 0.1 
ms (cr = 0.32 ± 0.02) for SGR 1806—20 events. The rgo emission times of bursts from both 

sources are on average shorter than Tgo durations. Note that the log-Gaussian mean values 
of Tgo distributions of single pulse bursts from both SGRs are quite similar to those of the 
Tgo distributions. 

3.3. Duty cycles (5go) 

As shown in the previous section, the rgo values of bursts are in general smaller than 
their Tgg values. The reason is that rgo is not a measure of the actual duration of the 
burst when the morphology of the burst is complex (e.g. multiple peaks). Mitrofanov et 
al. (1999) suggested that the ratio rgo/Tgo can then be used to describe a duty cycle, (5go. 
We, therefore, determined the 5go parameters for all bursts of both sources (Figure 7, left for 
SGR 1900+14 bursts, right ior SGR 1806-20 bursts). Both distributions peak at about 6go 
~ 0.45 and their overall shapes are very similar, indicating a strong similarity of the type of 
burst emission in the two sources, although their overall duration distribution differ. 

3.4. Duration - Fluence - Hcirdness correlations 

Gutenberg and Richter (1956) showed that there exists a power law relation between the 
magnitude (which is related to the total energy involved) of earthquakes and the durations 
of strong ground-shaking at short distances from an earthquake epicenter. We investigated 
whether a similar correlation exists in SGR events using their total burst counts and their 
Tgo durations (as estimated in §3.1), for both SGR 1806-20 and SGR 1900+14. For each 
set of bursts, we grouped the Tgo values into logarithmically spaced bins and determined the 
weighted mean value of total burst counts and Tgo durations for each bin. In Figure 8, we 
show the plot of integral counts vs Tgo durations of SGR 1806—20 bursts. The crosses in this 
figure are the weighted means of each parameter; the errors on the mean-Tgo values denote 
the range of each bin while the errors on the mean counts are due to sample variance. The 
dark points are the individual measurements of the burst counts and Tgo values of single 
pulse events and the gray circles are those of multi-peaked bursts. Figure 8 shows that the 
integral counts and durations of SGR 1806—20 bursts are well correlated. To quantify this 
correlation, we determined the Spearman's rank order correlation coefficient, p = 0.91, and 
the probability of getting this value from a random data set, P=3.4 x 10~"^. We further fit a 
power law model to the mean values (crosses) of the data using the least squares technique, 
which yields a power law index of 1.05 ± 0.16. 
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Similarly, Figure 9 shows that the integral counts and T90 durations of SGR 1900+14 
bursts are also correlated, having p = 0.89 and P = 8.6 x 10""^. A power law fit to the 
mean values of the data yields an index of 0.91 ± 0.07. The asterisk shown in the upper 
right portion of Figure 9 indicates where the precursor of the August 29 burst falls. This 
event was exceptionally long, bright and resembled the August 27 giant flare in various ways 
(Ibrahim et al. 2000). A spectral line at 6.4 keV was reported during the precursor of event 
of August 29 burst with ~ 4 cr significance (Strohmayer & Ibrahim 2000). If spectral hnes 
are event intensity dependent, there are very few events in Figure 9 during which we would 
expect to see any hnes. The detailed spectral analysis of all SGR bursts is underway. 

In conclusion we find a good correlation between the duration and total counts (energy) 
of SGR bursts, quite similar to the one estabhshed for earthquakes. It is noteworthy that 
the burst count fluences of single pulse events (from both systems) span a range almost as 
wide as that of fluences of multi-peaked bursts. 

In order to investigate burst spectral variations versus fluence, we calculated an event 
hardness ratio defined as the ratio of the total burst counts observed between 10—60 keV 
to those between 2 — 10 keV. Bursts with total counts less than ~ 50 yielded statistically 
insignificant hardness ratios (< 3 a), and therefore, were excluded from our analysis. For 

SGR 1806—20 bursts, wc have 159 events with hardness ratios measured to > 3 cr accuracy. 
We divided the total counts of these events into logarithmically spaced bins and determined 
the weighted mean hardness ratio for each group. Figure 10a shows that the SGR 1806—20 
burst hardness ratios arc anti-corrclatcd with fluence (p=— 0.96, P=2.6 x 10^'^). For 385 
SGR 1900+14 bursts. Figure lOd shows a marginal anti-correlation between hardness and 
fluence (p = —0.89, P=4.6 x 10~^). Interestingly, we flnd that although the hardness — 
fluence anti-correlation is evident in both sources, the SGR 1806—20 events are owera/Z harder 
than the SGR 1900+14 ones. Our energy selection optimizes the PGA energy response, so 
that we are not affected by instrumental biases. 

Fenimore, Laros & Ulmer (1994) have performed a similar study of hardness ratio vs 
fluence for 95 SGR 1806-20 events detected with IGE. We cannot directly compare our results 
with this study for two reasons. The energy ranges over which hardness ratios are computed 

differ significantly; while wc use (10 — 60 kcV)/(2 — 10 keV), Fenimore et al. (1994) use 
(43.2 - 77.5 keV)/(25.9 - 43.2 keV). The PGA sensitivity drops significantly beyond 25 
keV, limiting the possibility of comparison. Furthermore, our fluence range ends at 2 x 10^^ 
ergs cm~^, just below the range described in Figure 1 of Fenimore et al. (1994). However, 
since Fenimore et al. report a constant trend above 10~^ ergs cm~^, one could postulate that 
the hardness ratio vs fluence trend we report, levels off at higher fluences (energies). 

We further investigated the hardness — fluence trend for burst morphology sub-sets. 
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namely for single pulse and multi-peaked burst groups for each source. Figure 10b exhibits 
hardness — flucnce plot of single pulse SGR 1806—20 events and Figure 10c shows that of 
multi-peaked bursts. We see that both sets display spectral softening as the burst count 
fluence increases. Similarly SGR 1900-1-14 events are shown in Figure lOe and f. For this 
source, the hardness — fluence anti-correlation is signiflcant only for single pulse bursts. 
It is important to note that the peak rates of most of the highest fluence bursts reach > 
10^ counts s"^ (on 1/1024 s time scale) around which the PGA pulse pileup may become 
important. This effect can artificially harden the observed count spectrum at these rates. 
The spectral hardening or leveling off seen in the last bins of the plots of Figure 10 may well 
be due to the pulse pileup effect and should not be considered as intrinsic source property. 

We next investigated the relationship between the hardness and duration of SGR bursts. 
Although we find hardness — fiuence anti-correlations and fiuence — duration correlations, 
both SGR 1806-20 (Figure 11, squares) and SGR 1900+14 (Figure 11, diamonds) burst 
hardness ratios are independent of the event durations. 

4. Discussion 

Our study demonstrates that unhke the Tgo duration distribution of cosmic GRBs which 
shows a bi-modal trend with peaks at ~ 0.31 s and ~ 37 s (Kouveliotou et al. 1993; Paciesas 
et al. 1999), the Tgo distribution of SGR bursts displays a single peak which varies for each 
source: ~ 93 ms and ~ 162 ms for SGR 1900+14 and SGR 1806-20, respectively 

We find that the Tgo durations of single-pulse bursts from both SGRs form narrow 
distributions (compared to those of multi-peaked events) which peak at ~ 47 ms and ~ 88 
ms for SGR 1900+14 and SGR 1806-20, respectively. These bright, hard, and short bursts 
are almost certainly powered by a sudden disturbance of the rigid neutron star crust, which 
transmits energy to the magnetosphere. Magnetic stresses, which force the star between 
distinct metastable equilibria, provide the most plausible source of energy for the giant 
flares, and by inference for the short bursts (TD95; Thompson & Duncan 2001). 

The duration of the multi-peaked bursts is clearly flxed by the time between successive 
releases of energy. The widths of the single pulse bursts could, in principle, be limited either 
by the rate of release of energy from the initial reservoir or, alternatively, by the time for 
the released energy to be converted to radiation through some intermediate reservoir. There 
have been various suggestions for such an intermediate storage mechanism: a hot fireball 
that is confined on closed magnetic field lines (TD95); a region of strong magnetic shear and 
high current density (Thompson et al. 2000); or a persistent vibration of the star (Fatuzzo 
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& Melia 1994). 

Measurement of the burst rise time tr provides a discriminant between these possibihties. 
We find that tr is characteristically much shorter than the total duration, as defined by 
either T90 or T90. Moreover, the distribution of U/Tgo is broad (both for single-pulsed bursts 
and multi-peaked events), which suggests that tr is not directly connected to coohng. For 
example, the X-ray luminosity of a trapped fireball in local thermodynamic equilibrium is 
proportional to its surface area, and is a weak function of its internal temperature (Thompson 
& Duncan 1995). The characteristic radiative timescales are, in general, very short at the 
high spectral intensity of an SGR burst. Thus, the rise of the X-ray flux plausibly represents 
the initial injection of energy, although strictly tr only sets an upper bound to the injection 
timescale. 

While the overall shape of T90 (both single and multi-peaked) and the rgo distribu- 
tions of both sources are quite similar (relatively consistent Gaussian widths), the peaks 
of both distributions for SGR 1900-1-14 occur at shorter durations compared to those for 
SGR 1806—20. This systematic difference in the burst durations probably results from some 
differing intrinsic property of the sources, such as the strength of the magnetic field, or the 
size of the active region. 

We find a power law correlation between the total burst counts (fiuence) and the duration 
of SGR bursts with a power law index around 1. Similar behavior was noted for earthquakes 
by Gutenberg & Richter (1956). Within the context of earthquake mechanics, one way 
of defining duration (also known as "bracketed duration" [Bolt 1973]) is the time between 
the first and last 5% excesses of g (gravitational acceleration on Earth) by the threshold 
acceleration of strong ground shaking. Recently, Lay & Wallace (1995) presented the power 
law correlation between the seismic moment (oc energy) and duration of 122 earthquakes 
with energies between 3.5 x 10^"^ erg and 2.8 x 10^^ erg. The power law fit to these events 
yields an index of 3.03. 

An equally important constraint on the injection and cooling mechanisms comes from 
the anti-correlation between the hardness and fiuence of the SGR bursts. Although very 
significant for SGR 1806-20, this anti-correlation is much milder than that expected for 
black-body emission from a region of constant area; indeed, the trend of increasing hardness 
with lower fiuence is opposite to that expected for constant area emission. Two basic types of 
radiative mechanism could reproduce this trend. First, the emitting plasma could be in local 
thermodynamic equilibrium, which requires that its size (radiative area) should decrease at 
lower fluences. An alternative possibility is that the spectral intensity of the radiation field 
sits below that of a black body, and that the temperature of the emitting plasma is buffered 
within a narrow range. We consider each of these possibilities in turn. 
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An SGR burst can be parameterized by the rate of injection of energy into the mag- 
netosphere, Linj, and the volume V of the injection region. When Ljnj is large and V is 
small enough, it is not possible to maintain a steady balance between heating and radiative 
cooling. The deposited energy is locked onto closed magnetic field lines of the neutron star, 
in a "trapped fireball" composed of photons and electron-positron pairs (TD95). This kind 
of event will tend to have a soft spectrum, because the injected energy has thermalized, and 
the plasma remains in LTE very close to its photosphere. The rise time is comparable to the 
time over which energy is initially injected, but the decay is limited by the rate of cooling 
through a thin radiative surface layer, which contracts toward the center of the fireball. The 
declining fight curve of the 27 August 1998 giant flare can be accurately flt by such a model 
(Feroci et al. 2001; Thompson & Duncan 2001). 

A second burst from SGR 1900+14 on 29 August 1998 has been interpreted in this 
trapped flreball model (Ibrahim et al. 2000). The main 29 August burst had a much shorter 
duration than the giant fiarc (~ 3.5 s versus ~ 400 s). This bright component was followed 
by a much fainter pulsating tail (extending out to ~ 1000 s) which provides direct evidence 
for heating and compression of the neutron star surface by the fireball. In this model, the 
short duration and high luminosity of the bright component require that the trapped fireball 
had an approximately planar geometry, as would be expected if the energy were released 
along an extended fault (Ibrahim et al. 2000). 

As our results show, most single-pulsed bursts from SGR 1900+14 have a much shorter 
duration (40 times smaller) than the bright component of the August 29 burst. It remains 
unclear, therefore, whether the trapped fireball model also applies to these much more fre- 
quent events. If it does apply, then the narrowness of the Tgo distribution (compared with 
the wide range of measured fluences) also requires a planar geometry, because the cooling 
time is determined by the smallest dimension of the flreball. 

The radiative mechanism is somewhat different if the injection luminosity Linj lies below 
a critical value of ~ 10**^ km) ergs s~^ (assuming a spherical geometry; Thompson 

& Duncan 2001). When the compactness is that low, it is possible to maintain a steady 
balance between the heating of a corona of electron-positron pairs, and radiative diffusion 
out of the corona (which one deduces must be optically thick to scattering). The central 
temperature of the corona is buffered in the range ~ 20 — 40 keV, and remains higher at 
lower luminosities - just the trend observed for SGR 1806-20 (Fig. 10). If this radiative 
model applies to the majority of single-pulsed events, then the constraints on the geometry 
of the active region are weaker, and one must consider alternative mechanisms for storing 
the injected energy, such as a persistent current driven by shearing motions in the neutron 
star crust (Thompson et al. 2000). 
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Fig. 1. — {top panel) Light curve of a SGR 1900+14 burst (2—60 keV) accumulated with 20 
ms time resolution. The vertical dot-dash lines before and after the burst are the boundaries 
of the selected pre- and post-background intervals, {middle panel) The solid line shows the 
cumulative counts, with the background contribution (as fit by a first-order polynomial) 
subtracted off. The dashed fine shows a step-function fit to the burst counts. The lower 
horizontal dotted line represents the 5% and upper horizontal dotted shows the 95% count 
fiuence level. Note that the top and the middle panels are in the same time scale, {bottom 
panel) Time profile of the same burst with 7 ms time resolution and zoomed in near the time 
of the burst. The vertical dotted lines are the start and the end times of its Tgo duration. 
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Fig. 2.— Time profiles of some SGR 1900+14 and SGR 1806-20 bursts (on 7 ms time scale) 
to illustrate single pulse bursts (top plots) and multi-peaked bursts (middle and bottom 
plots) . 
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Fig. 3. — Distributions of the Tgo-durations for all (solid), single pulse (dashed) and multi- 
peaked (dash-dot) SGR 1900-1-14 bursts. The solid, dashed and dash-dot curves are obtained 
by fitting a log-Gaussian model to each set. Sets peak at 93.4 ms, 46.7 ms and 148.9 ms, 
respectively. 



10 



100 

Tgo (ms) 



1000 



Fig. 4. — Distributions of the Tgo-durations for all (solid), single pulse (dashed) and multi- 
peaked (dash-dot) SGR 1806—20 bursts. The solid, dashed and dash-dot curves are obtained 
by fitting a log-Gaussian model to each set. Sets peak at 161.8 ms, 88.1 ms and 229.9 ms, 
respectively. 
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Fig. 5. — Distributions of the ratios of rise times, t^ to Tgg durations for single pulse bursts 
(solid) and multi-peaked bursts (dashed) from SGR 1900+14 (left) and SGR 1806-20 (right). 
Note that the number of multi-peaked events with the lowest t^/Tgo ratios is artificially large 
due to the longer (compared to t^) durations of these bursts (see Figure 2, middle and bottom 
plots) . 
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Fig. 6.— Distributions of the rgo emission times for SGR 1900+14 (left) and SGR 1806-20 
(right) bursts. The dashed curves are obtained by fitting a log-Gaussian model to each set, 
with peaks at 49.6 ms and 82.3 ms, respectively. 
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Fig. 7.— Distributions of the duty cycles, Sgo, for SGR 1900+14 (left) and SGR 1806-20 
(right) bursts. The Gaussian fits peak at 0.45 and 0.46, respectively. 
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Fig. 8. — Plot of the total counts vs Tgo-durations for 190 SGR 1806— 20 bursts which shows 
a correlation (p = 0.91, P—3A x 10~^). The sohd line is obtained by fitting a power law 
model to the crosses (which are explained in the text). The best fit power law index is 1.05 
± 0.16. The circles are the individual measurements for single pulse (dark) and multi-peaked 
(gray) bursts. 
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Fig. 9. — Plot of the total counts vs the Tgo-durations for 455 SGR 1900+14 bursts which 
shows a correlation (p = 0.89 and P = 8.6 x 10~^). A power law fit to the crosses (solid 
line) yields an index of 0.91 ± 0.07. The asterisk indicates where August 29 precursor would 
be seen. The circles represent the measured values for single pulse (dark) and multi-peaked 
(gray) bursts. 
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Fig. 10. — Plot of hardness ratios vs total counts for (a) all, (b) single pulse, (c) multi- 
peaked SGR 1806—20 bursts. Similarly, those for (d) all, (e) single pulse, (f) multi-peaked 
SGR 1900+14 events. 
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Fig. 11. — Plot of hardness ratios vs Tgo durations for SGR 1806—20 bursts (squares) and 
for SGR 1900+14 (diamonds). 



